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ABSTRACT 

Aims. The Large Binocular Cameras (LBC) are two twin wide field cameras (FOV ~ 23' x 25') mounted at the prime foci of the 8.4m 
Large Binocular Telescope (LBT). We performed a weak lensing analysis of the z = 0.288 cluster Abell 61 1 on g-band data obtained 
by the blue-optimized Large Binocular Camera in order to estimate the cluster mass. 

Methods. Due to the complexity of the PSF of LBC, we decided to use two different approaches, KSB and Shapelets, to measure the 
shape of background galaxies and to derive the shear signal produced by the cluster. Then we estimated the cluster mass with both 
aperture densitometry and parametric model fits. 

Results. The combination of the large aperture of the telescope and the wide field of view allowed us to map a region well beyond 
the expected virial radius of the cluster and to get a high surface density of background galaxies (23 galaxies/arcmin^). This made 
possible to estimate an accurate mass for Abell 611. We find that the mass within 1.5 Mpc is: (8 ± 3) x 10'"* Mq from the aperture mass 
technique and (5 ± 1) x 10'*Mg using the model fitting by a NFW mass density profile, for both Shapelets and KSB methods. This 
analysis demonstrates that LBC is a powerful instrument for weak gravitational lensing studies. 

Key words. Cosmology: dark matter; Galaxies: clusters: individual; Gravitational lensing 

' 1. Introduction estimates. 

Several techniques are commonly used to estimate masses for 

According to the hierarchical model of structure formation, clus- galaxy clusters: the X-ray luminosity or temperature of the hot 

ters of galaxies are the most massive objects in the universe mtracluster gas, the Sunyaev-Zel'dovich effect, the number of 

and the cluster mass function is a powerful probe of cosmolog- bright galaxies in a cluster, and the velocity dispersion of the 

ical parameters (e.g. 'Evrard' 1989': Eke et al■ ^l 9981: lHenrvll2000l: cluster galaxies. The disadvantage of all of these methods is that 

[Allen et al. 2004; Vikhlinin et al. 2009). In addition, the ratio be- *ey are mdirect and requn-e significant assumptions about the 

tween the cluster gas mass, as estimated with X-ray observa- dynamical state of the cluster Gravitational lensing, in conti-ast, 

tions, and the total mass in a galaxy cluster provides stringent is only sensitive to the amount of mass along the Une of sight and 

constraints on the total matter density. Specifically, the apparent allows reconstruction of the projected cluster mass regardless 

evolution of the gas fraction with redshift can be used to esti- its composition or dynamical behavior (e.g. [Kaiser & Squires | 

mate the contri bution of th e dark energy component to the cos- 1993). The only direct method to estimate cluster masses is 

mic density (e.g lAUen et al.l l2008: Ettori et al. 2009). The use of therefore via measurement of the distortion (shear) of the shapes 

clusters as a cosmological probe therefore requires reliable mass of background galaxies that are weakly lensed by the gravita- 
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Fig. 1. A three-color image (2' x 1.6') of Abell 611 obtained 
with LBC observations in the u-, g-, and r-bands. A giant arc 
is clearly visible close the brightest cluster galaxy {a= 08h 00m 
57s,5=+36d 03' 23"). 



tional potential of the cluster. 

This distortion is very small and lensed galaxies are usually at 
high redshift. Observational studies to measure weak gravita- 
tional lensing by clusters require deep images in order to detect 
these faint sources and to obtain a high number density of back- 
ground galaxies. Moreover this kind of analysis requires very 
high quality images in order to measure the shape of the lensed 
sources with high precision: good seeing (< 1") conditions and a 
high signal-to-noise ratio (typically > 10) are needed. Wide-field 
images are also required to obtain a statistical measure of the tan- 
gential shear as function of distance from the cluster center so 
that the projected mass measured by weak lensing includes es- 
sentia lly all of the mass of the cluster (IClowe & Schneideil200ll 
I2002h . 

In the last decade substantial progress has been made with weak 
lensing studies thanks to the advent of wide-field data with linear 
detectors, the development of sophisticate d algorithms for shape 
measurer nents ( e.g. Kaiser et al. 1995, Be rnstein & Jarvisll2002l 
iRefr egier 2003, Kuiiken 2006), and the availability of multi- 
band photometry which provides information abou t the redshift 
distribution of the lensed sources (lllbertet al.ll2"006l) . 
Here we describe the results of a weak lensing analysis of the 
Abell 61 1 cluster This analysis is based on images obtained with 
the Large Binocular Camera (LBC), which are a pair of prime 
focus cameras mounted on the two 8.4m diameter mirrors of the 
Large Binocular Telescope (LBT). Each LBC has a 23' x 25' 
field of view (FOV) and, combined with the collecting area of 
LBT, is a very powerful instrument for weak lensing studies. 
Abell 611 is a rich cluster at redshift z = 0.288 
(IStruble&Herbertll999h that appears relaxed in X-ray data, has 
a regular morphology, and the brightest cluster galaxy (BCG) 
is coincident with the center of X-ray emission (Donnarumma 
et al. in preparation). A giant arc due to strong lensing is also 
clearly visible close to the BCG (Fig.[TJ. 

In this work we describe a weak lensing analysis to estimate the 
mass of Abell 611 from a deep g-band LBC image whose field 
of view extends well beyond the expected virial radius of the 
cluster We compare the mass estimated from gravitational lens- 
ing with previous lensing results and with other mass estimates 
available in the literature that were derived by secondary tech- 



niques. In particular we compared our weak lensing results with 
new mass measurements obtained by X-ray analysis of Chandra 
data provided to us by Donnarumma et al. (in preparation). 
Mass measurements by weak lensing do not need any assump- 
tion about the geometry of the cluster; however, assumptions are 
required to compare projected lensing masses with other mass 
estimates. Thus projection effects have to be taken into account 
durin g this kind of analysis, in particular the true triaxiality of the 
halo foe Filippis et al.ll2005L lG avazz 1 12001 and the presence o f 



unrelated structures along the line of sight (lMetzleretal.ll200lL 
Hoekstra 2007) can be sources of noise or bias on the projected 
mass measurements. 

The paper is organized as follows. In the first sections we 
describe the data (Section |2]i used for a weak lensing anal- 
ysis of Abell 611, the catalog extraction of the background 
sources (Section[3]l, and the selection of candidate cluster galax- 
ies (Section nil. The two different pipelines used to extract the 
shear signal from the images are described in Section|5] and their 
results are compared in Section 15.31 Finally, both shear maps 
are used to estimate the mass of the cluster with different tech- 
niques (Section|6]l. The results are summarized and discussed in 
the Section |7] 

Throughout this paperwe adopt //(>= 70 km Mpc"', Q„,=0.3, 
and Qa=0.7. At the distance of Abell 611, 11' radius corre- 
sponds to a projected physical distance of nearly 3 Mpc. 

2. Observations and Data Reduction 

Abell 611 was observed in March 2007, during the Science 
Demonstration Time (SDT) for the blue-optimized Large 
Binocular Camera (LBC), which is one of the two LBCs built 
for the prime foci of the LBT. The LBC focal plane consists 
of four CCDs (2048 x 4608 pixels, pixel scale 13.5 fim, gain 
~ 2 e /ADU, readout-noise ~ 11 e ). The CCDs are arranged 
so that three of the chips are butted along their long edges and 
the fourth chip is rotated counterclockwise by 90 degrees and 
centered along the top of the other CCDs (see Fig|2]i. The field 
of view is equivalent to 23' x 25' and provides images with a 
sampling of 0.225"/pixel. Because each LBC is mounted on a 
swing arm over the primary mirror, the support structure lacks 
the symmetry of most prime focus instruments. Moreover LBC 
PSFs are dominated by optical aberrations from mis-alignments, 
which can cause PSFs to not have bilateral symmetry. This is 
a potential complication for the weak lensing analysis and we 
discuss this point further be low. More detail s about the charac- 
teristics of LBC are given in lGiallongo et al.l (120 08). 
The observations, collected in optimal seeing conditions 
(FWHM ~ 0.6"), consisted of several sets of exposures of 5 
minutes each in a wide M-band, SDSS g- and r-band filters. The 
total exposure time was 1 hour in g, 15 minutes in r, and 30 min- 
utes in M. Unfortunally some M-band observations were not us- 
able, so the real total exposure time used for this band is 20 min- 
utes. For the present work we used the deep, g-band data for the 
weak-lensing analysis and the u- and r-band data to select cluster 
galaxies. Each image was dithered by 5 arcseconds in order to 
remove bad pixels, rows, columns, and satellite tracks. This off- 
set is not large enough to fill the gaps between the CCDs, but the 
analysis plan was to treat CCDs separately due to expected PSF 
discontinuities at the chip boundaries. The offsets were therefore 
kept small to maximize depth and uniformity. The images were 
reduced by the LBT pipelin43 implemented at INAF-OAR. The 
flat-field correction was done using both a twilight flat-field and 
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Fig. 2. g-hand image of the full field of LBC, centered on Abell 
611. The box marks the region of 5000 pixels we used for the 
analysis. 
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Table 1. Exposure times, limiting magnitudes for point-like 
sources and zero points (AB) for the observations in each band. 



a superflat obtained during the night. Moreover a geometric dis- 
tortion correction was performed in order to normalize the pixel 
size, which showed differences across the CCDs due to field dis- 
tortions in the optics. The astrometri c solution was computed 
using the ASTROMC package (Rad ovich et alj|2008l) . This so- 
lution was then used to resample and coadd the images using the 
SWariS software. 

Standard fields for photometric calibration were not observed 
during the SDT. So we used the values of zero points for each 
band (Tab. [B given by iGiallongo et"an (l2008h . Table □ also 
shows the limiting magnitudes estimated from the faintest point- 
like objects detected at the 5cr and 20cr level. 

3. Catalog extraction 

The dete ction of sources was pe rformed using the SExtractor 
package (iBertin & ArnoutsI 1996b . Regions of the image present- 
ing potential problems, such as spikes and halos around bright 
stars, were masked by visual inspection; sources inside such re- 
gions were discarded from the final catalog. In addition, we re- 
moved sources located at the borders of each CCD, where the 
SNR was lower due to the small dither offset. Finally, a very 
bright star dominates one of the CCDs. We therefore decided to 
limit our analysis to a box (displayed in Fig|2]i with a size of 
5000 pixels (corresponding to ~ 18.7') centered on Abell 611. 

^ developed by E. Bertin, [http://terapix.iap.fr| 
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Fig. 3. Magnitude (g) vs. half-light radius (r;,) plane. Green zone 
shows the background galaxies selected for the lensing analysis 
(23 < g < 26, rh > 1.8 pixels) and pink zone the unsaturated 
stars selected for correction of PSF anisotropy (20 < g < 23, 
1.4 < r/, < 1.8 pixels). See |3]for further details. 



Starting from this box of 350 arcmin^, the eff'ective area used for 
the analysis was ~ 290 arcmin^ after removing all the masked 
regions (30% due to bright stars, 70% due to regions between 
adjacent CCDs with no data or low S/N). 
The separation between stars and galaxies was performed in the 
mag - r/, plane, where magnitudes (mag) and half-light radii (r/,) 
were obtained from the MAG JVUTO and FLUX_RADIUS pa- 
rameters computed by SExtractor. Unsaturated stars were se- 
lected on the vertical branch (see Fig. [3]) in the range 20 < ^ < 23 
mag and 1.4 < r/, < 1.8 pixels. In this way we obtained 302 stars 
for the PSF correction, with a SNR> 200 for the faintest ones. 
For the lensing analysis, only background galaxies located at 
redshifts larger than z - 0.288, the redshift of the cluster, should 
be used. Unfortunately, in our case the number of available bands 
does not allow us to estimate accurate photometric redshifts of 
these faint galaxies. The selection of the background galaxies 
was therefore done by choosing an adequate cut in apparent 
magnitude. 

The choice of the upper magnitude limit was ba sed on the galaxy 
redshift distribution obtained by lllbert et al] 1*2006) from the 
Canada-France-Hawaii Telescope Legacy Survey (CFHTLS), 
which also used the SDSS photometric system. Taking into 
account th e accu racy of the photometric redshifts (3%) of 
Ilbe rt et all (|2006), the approximations due to the different bands 
they used compared to ours and the assumption that we have the 
same galaxy distribution in our field, we chose 0.4 as redshift ref- 
erence value to perform the magnitude cut. Galaxies with z < 0.4 
were assumed to belong to the cluster or be foreground galaxies. 
This reference value was chosen to be larger than cluster red- 
shift, in order to reduce the contamination of foreground galax- 
ies as much as possible, taking into account the approximations 
discussed above. 

Figure |4] shows the fraction of the total CFHTLS sources at 
z < 0.4 (upper panel) and the fraction of background galaxies 
at z > 0.4 (bottom panel) as a function of the apparent mag- 
nitude cut. From this figure we conclude that a magnitude cut 
at g > 23 is a good compromise to minimize the contamination 
from likely foreground and cluster galaxies (~ 10%) and to max- 
imize the number density of background galaxies (~ 98%). The 
faint magnitude cut was chosen at ^ < 26, that is the magnitude 
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Fig. 4. The fraction of likely foreground and cluster galaxies 
(top-panel) and background galaxies (bottom-panel) from the 
CFHTLS versus magnitude cut in g-band. The dotted vertical 
line is the lower limit of the magnitude cut adopted to select 
background galaxies. See |3]for details. 



limit where we have a signal-to-noise ratio SNR > 10 for the 
sources, where SNR is defined as FLUX/FLUX_ERR as mea- 
sured by SExTRACTOR. The final catalog contained 8134 back- 
ground galaxies. 

4. Candidate cluster members 

The candidate cluster members were selected from the simul- 
taneous usage of U-, g- and r- band photometry. To this end 
catalogs were extracted from these images, running SExtractor 
in dual-mode with the g-band image as detection image. Only 
sources detected in all three bands were used. 

We then applied an algorithm ( Fu et al., in prepara tion) simi- 
lar to the C4 Clustering Algorithm (Miller et al. 2005). The algo- 
rithm is based on the assumption that galaxies in a cluster should 
have similar colors and locate together in space. It evaluates the 
probability of each galaxy to be field-like: hence candidate clus- 
ter galaxies are those for which this probability is below a given 
threshold, as outlined below. 

1. Each galaxy was set in a four-dimensional space of a, 6, 
u - g and g - r. For each galaxy (named "target galaxy"), 
we counted the number of neighbors within the four- 
dimensional box, Maiget- The angular size of the box was set 
to 1 h'^ Mpc (~ 5.5' at the redshift of Abell 611). The sizes 
of the boxes in two color dimensions were determined as 



(5m« = ycr2^(stat)-Ho-2,,(sys), (1) 

where cr„,„(stat) is the observed error for two magnitudes 
(m, n), and cr,„„(sys) is the intrinsic scatter of the color m-n. 
For LBC A61 1 data, dug, 5gr are 0.49 and 0.31 respectively. 
2. This four-dimensional box was placed on 100 randomly 
chosen galaxy positions and at each position we counted 
the number of neighbors. These randomized number counts 
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Fig. 5. Colour- magnitude plot of the galaxies in the Abell 611 
field. Red points are the candidate cluster galaxies selected by 
the C4 method. Also displayed (solid line) is the results of the 
biweight regression fit. 

constructed a distribution of counts for the target four- 
dimensional box. This distribution is represented by the me- 
dian value A^median of the randomization counts. 

3. The probability p of the target galaxy being field-like was 
derived by comparing the target galaxy count A^target to the 
distribution of randomization values A^median- 

4. The distribution of p values was derived by repeating the 
above steps for all galaxies. We ranked the p values from 
smallest to highest and derived the value after which p starts 
to rise significantly. In this way we identified ~ 150 galaxies 
at r < 23 mag as the candidate cluster members. We further 
removed outliers in the g - r \s u - g diagram, leaving 125 
candidate members. 

Starting from these galaxies, we fitted the g - r vs. r Red 
Sequence (g - r = a -\- b ■ r) using a biweight regression method 
(Fig.|5]i, and obtained; a = 2.39+;;;^', b = -0.04+j; 

5. Weak lensing analysis 

Weak lensing is based on the measurement of the coherent dis- 
tortion of the shapes of background galaxies produced by a dis- 
tribution of matter This distortion is very small and it is un- 
measurable for any single background galaxy because the galaxy 
intrinsic ellipticity is not known. The dispersion of intrinsic el- 
lipticity is a source of noise which is oc cr^/ -^V- A statistical 
approach is therefore required, where the distortion can be mea- 
sured for a large number of sources in order to bring down that 
noise. This requires a careful treatment of systematic effects, as 
the shapes of galaxies may also be affected by contributions to 
the point spread function (PSF) by both the telescope and the 
atmosphere. In the last decade several methods have been devel- 
oped for this kind of analysis. The most popular is the KSB ap- 
pro achj_originaUy_prOTOsed^ et al. ( 1995) and improved 
by Luppino & Raised (Il997h and lHoekstra et al., (.1998.) . Several 
different implementations of this method exist in the literature 
and have been used in a number of KSB analysis pipelines. 
More recently, Refregier (2003) and lMasse y & Refregier (200j) 
have proposed a new method based on Shapelets. Several avail- 
able pipelines have also used this approach to measure the shear 
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signa l in various ways (e.g. lMassev & Refregieill2005L iKuiikenI 
I2006h . 

One of the key differences between these two approaches is the 
treatment of the PSF. The KSB method assumes that the PSF 
can be written as a convolution of a very compact anisotropic 
kernel with a more extended, circular function. These two terms 
are expressed in terms of the quadrupole moments of the sur- 
face brightness. As summarized below, the isotropic component 
is subtracted from the measured ellipticity and the anisotropic 
component is subtracted from a responsivity term. In contrast, 
in the Shapelets approach there is no assumption about the PSF 
shape. Individual galaxy images are decomposed into a com- 
plete orthonormal basis set consisting of Hermite (or Laguerre) 
Polynomials and the PSF correction is performed through de- 
convolution. 

Since the PSF of LBC presents a significant deviation from sym- 
metry, the analysis of Abell 611 presents a very good opportu- 
nity to compare the results produced by these two methods. For 
this comparison we started with the same initial catalogs of stars 
and galaxies for both pipelines. Specifically, we only considered 
sources with a SExtractor FLAG < 4, which removes sources 
that are possibly blended. As the subsequent steps performed by 
each algorithm are different, the same galaxy may be rejected 
by one algorithm but not by the other. The result is that dif- 
ferent output catalogs are produced by the KSB and Shapelets 
pipelines. To have a homogeneous comparison, we therefore fi- 
nally selected only the sources common to both output catalogs. 
Sources with an unphysical ellipticity \e\ > 1 were also removed 
from the final catalogs. 

5.1. KSB method 



We u sed the weak lensing pipeline described in iRadovich et al.l 
(12008) to compute the quantities relevant to the lensing analysis. 
This pipeline implements the KSB approach using a modified 
version of the IMC AT PI tools that was provided t o us by T. Erben 
dErben et al.ll200ll and lHetterscheidt et al.ll2007l) . 
In the KSB approach stars and galaxies are parametrized accord- 
ing to the weighted quadrupole moments of the intensity distri- 
bution using a Gaussian weight function whose sc ale length is 
the si ze of the source (the formalism is described in lKaiser et al.l 
Il995h . The main assumption of this approach is that the PSF can 
be described as the sum of a large isotropic component (seeing) 
and small anisotropic part. In this way the observed ellipticity 
Cobs can be related to the intrinsic source ellipticity e, and shear 
7 by the relation: 



^ohs 



(2) 



where f ™' is the smear polarizability tensor and f ^ is the pre- 
seeing polarizability tensor, which is related to t he shear polar- 
izabi lity tensor f and f ™ by Eq. |5] (see also iHoekstra et al.l 
11998 ). The quantity p characterizes the anisotropy of the PSF 
and is estimated from stars, which have zero intrinsic ellipticity: 



P = 



obs 



(3) 



If we average over a large number of sources, assuming a ran- 
dom orientation of the unlensed galaxies, we expect (e ,) = and 
so 

r - <^>, (4) 
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Fig. 7. P '''* / P""* values computed in different bins of r/,, for each 
CCD. These values are lower for the central CCD (#2) where the 
quality of PSF is better. 



where e,jo = Cobs - P' 
anisotropic distortions. 



is the ellipticity corrected for 



We computed the contribution due to the PSF anisotropy (Eq. 
|3]l from the selected stars. This quantity changes with position 
in the image, so we needed to fit it on each CCD in order to 
extrapolate its value at the position of the galaxy we want to 
correct. In our case we performed this fitting on each CCD and 
a second-order polynomial fit was sufficient. 
Fig.|6]displays the spatial pattern of the ellipticities for the stars 
in all CCDs, before and after the PSF correction. 
After that we computed P^ for each source: 



P^ = P 



sh 



psh* 



(5) 



http://www.ifa.hawai.edu/~kaiser/imcat 



As the PSF correction was done separately on each CCD, we 
decided not to fit as a function of the coordinates and instead 
we took an average value. 

Stellar ellipticities should be co mputed using the s ame 
weight function used for galaxies (see iHoekstra et al.l[T998h . so 
we considered a sequence of bins in r/, and for each galaxy we 
selected the PSF correction terms computed in the closest r;, bin. 
Fig.Qshows the values computed in different bins of r/,, for 
each CCD. 

We weighted the shear contribution from each galaxy according 
to: 

w = (6) 

as inlHoekst ra et al.l (12000*). where <Teo ~ 0.3 is the intrinsic rms 
of galaxy ellipticities, (Ae-)'''^ is the uncertainty in the measured 
ellipticity, where the average is computed on bo th components 
of ellip ticity for each galaxy (see eq. A8, A9 in IHoekstra et al.l 
(I2OOOI) '). 

A crucial point in this kind of study is the selection of the galax- 
ies to use for the shear analysis, as the contamination of fore- 
ground galaxies can dilute the lensing signal and lead to an un- 
derestimate of the mass. 
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Fig. 6. Removal of PSF anisotropy in the KSB approach. The observed (top-left), fitted (top-right) and residual (bottom-left) ellip- 
ticities of stars for all CCDs. The observed (black points) versus corrected (green points) ellipticities are also shown (bottom-right). 



Since the PSF degrades somewhat at the borders of the image, 
we limited our analysis to an 18.7' x 18.7' box centered around 
the BCG, as discussed in Section [3] After that we filtered the 
source catalog using the following criteria: P'' > 0.25, SNR 
> 10, r/, > 1.8 pixels, 23 < g < 26, ellipticities smaller than 
one, obtaining a surface density of ~ 25 galaxies/arcmin^. 
The cut P'' > 0.25 allowed us to discard sources that appeared 
too circular (e.g. stars incorrectly classified as galaxies). We con- 
sidered only galaxies with SNR > 10 in order to avoid noisy ob- 
jects, which can be a source of error in the computation of the 
shear signal. Finally, we used the magnitude cut 23 < g < 26 
to select background galaxies, whose choice was previously ex- 
plained in Section [3] 

5.2. Shapelets 

Another approach for weak lensing analysis is the use of 
Shapelets, which are basis functions constructed from two- 
dimensional Hermite polynomials weighted by a Gaussian. The 
translation, magnification, rotation and shear of astronomical 
images can be expressed as matrices acting on Shapelets coeffi- 
cients. The advantage of Shapelets is that a galaxy image can be 
described in reverse order: pixelation, convolution with the PSF, 
and finally distortion by shear. Shapelets have a free scale ra- 
dius j6 which is the size of the Gaussian core of the functions. Its 
truncated expansion describes deviations from a Gaussian over 
a particular range of spatial scales, which widens with order A^. 
For the analysis of the LBC data of Abell 611 we used the 



Shapelets pipeline developed by iKuiikenI ( 12006 ^. starting with 
the same stars and galaxies used in the KSB method. Each star 
was first fit with a circular Gaussian and the median radius was 
computed. This radius was then multipl ied by a f actor o f 1 .3 for 
the Shapelets fits, which was found bv IKuiikenI (l2006h to work 
well for a range of model PSFs up to Shapelets order N - ^. 
Then we obtained a Shapelets description of the PSF for each 
star In order to estimate the PSF model at the position of each 
galaxy, the Shapelets coefficients were interpolated by a fourth- 
order polynomial on the whole image frame (Fig.[8]l. The resid- 
ual of the PSF model fitting is shown in Fig.|9l 
The ellipticity of each source is then determined by least-squares 
fitting a model, which is expressed as the shear applied to a cir- 
cular source to fit the object optimally. The extension order of 
Shapelets for galaxies is taken the same as for stars. Performing 
the least-squares fit, the minimu m of can be fo und in a few 
Levenberg-Marquardt iterations dPress et al.lll986l) . The errors 
of Shapelets coefficients for each source are derived from the 
photon noise, and these can be propagated through in the 
function. Thus the error of shear measurement o-y is calculated 
from the covariance matrix which is given by the second partial 
derivatives of^^ at the best fit. 

The shear contribution from each galaxy is weighted according 
to: 



which combines the error in the shear measurements o-y and the 
intrinsic scatter cr^o- The last quantity was computed from the 
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X= 602 1130 2845 3915 5077 6332 

Fig. 8. Shapelets PSF models interpolated at different positions on each CCD. The distribution of these models corresponds to their 
actual placement on the CCD mosaic, CCDl to CCD3 from left to right in the bottom panel and CCD4 in the top. Contours are the 
representations of PSF shape decomposed by Shapelets. The X and Y values correspond to the pixel position. 



distribution of ellipticity components of galaxies. 

Finally, we removed galaxies with SNR < 10 and that failed the 

Shape lets expansion and radial profile cuts defined by Kuiike^ 

(l2006l) . This eliminated galaxies that were not well detected or 

measured, providing a number density of ~ 26 galaxies/arcmin^. 

Further details about the selection criteria are available in 

lKuiikenl(l2006l) . 



5.3. Comparison between KSB and Shapelets ellipticities 

After matching the two output catalogs, our final background 
galaxy sample has a surface density of ~ 23 galaxies/arcmin^. 
In weak lensing studies, the background number density 
of ground-based te lescope usually is around 15 to 20 
galaxies/arcmin^ (e.g. lPaulin-Henriksson et al.l2007l) . It depends 
not only on the size of telescope, the exposure time, the color fil- 
ter, seeing condition, but also on galaxy selection criteria. 
The maps of the PSF correction computed by both methods show 
a good-quality PSF in the central regions that then degrades fur- 
ther from the center of the field. Nevertheless the final correction 
is good with fluctuations in the PSF anisotropy less than 1%. 

Figure[TO]compares the first component of the shear ji mea- 
sured by the KSB and Shapelets methods. It shows good agree- 
ment in the range of -0.5 < yi < 0.5. Some scatter is present 



0.5 



-0.5 



-0.5 





71 (KSB) 



0.5 



Fig. 10. Comparison between the first component of the shear 
from the Shapelets and KSB methods (results with the second 
components are very similar). 
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for very elongated galaxies, but the fraction of these galaxies is 
less than 5% in the final, common catalog. As shown in Fig.fTTl 
these strongly elongated galaxies are down-weighted nearly by 
a factor of 2 compared to small ellipticity galaxies, so that they 
do not affect the final mass measurements of the cluster (as de- 
scribed further below). A similar behavior is seen for the second 
component 72- These matched catalogs were used to compute 
estimates of the cluster mass as described in the next sections. 

6. Mass Measurements 

The relationship between the shear y and the surface mass den- 
sity is: 

7(6*) = - r d^0'D{e-9')K{ff), (8) 

where: D - jg^^g^^, k = 'Z/I.cr is the convergence, I.cr - 

4^ 5^ ~ "Sg 'Bfi critical density of the cluster, and D j, D/ 
and Dis are source-observer, lens-observer and lens-source dis- 
tances respectively. While weak lensing actually measures the 
reduced shear g - for k <^ \, g ~ y. 

An optimal approach for computing the mass requires knowl- 
edge of the redshift of each background galaxy. As we do not 
have this information, we assumed that the background sources 
all lie at the same redshift according to the single sheet approx- 
imation (King & Schneider 2001). An estimate of the redshift 
value to use for the weak lensing analysis was computed us- 
ing the first release of photometric redshifts available for the Dl 



1.5 1— 
1.4 -J 



0.8 - i „ - 

0.7 - ^ T - 

0.6 - I - 

0.5 I ' ' ' ' ' ' ' '-^ 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

|eil 

Fig. 11. The plot of the normalized galaxy weight as function 
of the absolute first component of ellipticity for KSB measure- 
ment. The applied normalization was Sw, x (5e 1 - 1 , where 5e 1 is 
the constant bin width of ei. The Shapelets measurement shows 
similar down-weight for elongated ellipticity galaxies. 



deep field of the CFHTLS, adopting the magnitude cut in the ^- 
band chosen here for the background galaxies selection (23 < 
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g < 26 mag; see Section[3]l and assuming a Gamma probability 
distribution (iGavazzi et al.l2004l) . This yielded a median redshift 
z - 1 .05 and for our analysis we assume all of the background 
galaxies lie at z ~ 1, correspondin g to (13) ={Di JDs) ~ 0.65. 
As discussed by iHoekstral (|2007|) . the single-sheet approxima- 
tion results in an overestimate of the shear by a factor 



1 -I- 



- 1 



(9) 



We computed this factor using the CFHTLS catalog of photo- 
metric redshifts and obtained = 1.167. 

The convergence k gives an estimate of the surface mass 
density apart from an unknown additive constant - the so-called 
mass-sheet degeneracy. We tried to solve this degeneracy using 
two different approaches: assuming either that k vanishes at the 
borders of the image or a particular mass profile whose expected 
shear profile is known. 

6.1. S-Maps 

In Figure[T2]we plot the so-called S-maps (ISchirmer et al.ll20(j4l) 
for these data. S-maps are computed as the ratio S = Map /cr Map 
where: 



M, 



ap 



^ie,,iWiQ{\9i - gpl) 



■Oo\) 



2(E,-w,)2 



(10) 
(11) 



For this calculation the image is considered as a grid of points, 
e, , are the tangential components of the ellipticities of the lensed 
galaxies, which are computed by considering the center of each 
point of the grid, w, is the weight as defined in Equation|6l and Q 
is a window function, chosen to be a Gaussian function defined 
by: 

Qi\o-eo\) = -^gxpf- ^^~f"^ ] (12) 



0? 



where Oq and 0c are the center and the size of the aperture. The 
ratio S - Map/cTM^p provides an estimate of the SNR ratio of 
the dark matter halo d etection. S-maps are discussed further by 
ISchirmer et an(l2004l) . 

We computed these maps using shear catalogs obtained from 
both the KSB and Shapelets pipelines. In both cases (see Fig. 
\T% the maps show that the lensing signal is peaked around the 
BCG, confirming that this is indeed the center of the mass distri- 
bution. The mass distribution also appears quite regular, which 
is in agreement with what is indicated by the X-ray maps. In 
Figure[T3]S-map contours are overlaid on the r-band luminosity- 
weighted density distribution of the red sequence galaxies of 
Abell 611, selected in Section |4] showing that the mass distri- 
bution follows that of the red cluster galaxies. 



6.2. Aperture densitometry 

In order to trace the ma ss profile of the clust er, w e computed the 
^-stati stic described in IClowe et"an (Il998h and iFahlman et al.l 
n99l : 

m )^-K(6<6i) - m <0< 0,„a,) = 2 {yT}d ln0 (13) 

Jfl, 

2 




{7T)dln0, 



Fig. 13. r-band luminosity-weighted density distribution of red 
sequence galaxies of Abell 611. The overlaid contour (black 
lines) is the S-map (computed and discussed in Section 16.11 ) 
showing the SNR of the shear signal around the cluster ob- 
tained from the Shapelets analysis. The levels are plotted be- 
tween o-jnin = 3.5 and o-^ax - 5. 



Mass measurements are computed within different apertures of 
increasing radius using a control annulus far from the center of 
the distribution (the BCG). Map - 7T0^{(0i)l,cr is the mass within 
the last aperture and provides a lower limit on the mass, unless 
the value in the control annulus is equal to zero. This method 
allows us to choose the size of the annulus that satisfies the de- 
sired condition; moreover it has the advantage that this mass 
computation does not depend on the mass profile of the clus- 
ter (Clowe et al. 1998). 

We chose 30" < < 500" and 0,^_, = 600", which yielded 
projected masses within ~ 1500 kpc of 7.7 ± 3.3 x IO^^Mq and 
8.4±3.8x lO'^^Mo using the KSB and Shapelets shear catalogs, 
respectively. 

Unfortunately we could not extend our analysis further from the 
center of the cluster because of the presence of a very bright star 
in the field that made the outer regions unusable. 
Since weak lensing is sensitive to the total mass along the line 
of sight, the observed aperture mass is the sum of the mass of 
the cluster and any contribution from other, uncorrected struc- 
tures along the line of sight. This contribution is assumed to be 
negligible in the central regions of the cluster, which are much 
denser, an d become more relevant in its outer regions. As dis- 
cussed by IHoekstral (12001). the effect of this contribution does 
not introduce any bias, but it does add a source of noise to the 
lensing mass. Aperture densitometry is more affected by this un- 
certainty than parametric methods because it is sensitive to the 
lensing signal at large radii. Nevertheless, for observations of 
rich clusters at intermediate redshifts, this uncertainty is fairly 
small because the bulk of the background sources are at much 
higher redshifts than the cluster. 

6.3. Model fitting 

The model fitting approach for estimating the mass of a cluster 
consists of assuming a particular analytic mass density profile 
for which calculate the expected shear and then fitting the ob- 
served shear with the model by minimizing the log-likelihood 
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120.26 120.24 120.22 120.26 120.24 120.22 

Fig. 12. S-maps obtained from the KSB (left panel) and the Shapelets (right panel) analysis. The levels are plotted between cr, 
3.5 and osnux - 5. They are overplotted on a g-band grey scale image (~ 4 arcmin) of the center of the field of Abell 611. 



function ([Schneider et al.ll2000l) : 



h)\' 



crHgiOi)] 



+ 2ln cr[gm 



(14) 



with cr[gm^ (I -g(ed^)cr,. 



In this analysis we assumed both a Singular 
Isothermal Sphere (SIS) and a Navarro-Frenk- White 
dNavarro. Frenk & Whitdl 19971) model. 

In the SIS model the density profile depends on one parameter, 
the velocity dispersion cr; 



Ps!s(r) = 



cr 



1 



27tG r2 ■ 

In this profile the shear is found to be related to cr by: 

In (T^ Du Oe 



jriO) 



c2 D, 



(15) 



(16) 



(e.g. lBartelmann & Schneider 

The mass density profile predicted by the Navarro-Frenk- White 
model (hereafter NFW) is: 



PNFw(r) - 



ScPc 



(r/r,)(l+r/r,)2 



(17) 



where pc - 3H^{z)/i?inG) is the critical density of the universe at 
the cluster redshift, rs is a scale radius related to the virial radius 
by means of the concentration parameter Ci,,> - ryirlrs and 6c is 
a characteristic overdensity of the halo: 



6c 



3 ln(l -he)- c/(l H-c) 



(18) 



where A,,,> is the virial overdensity, approximated by A,,,,- ~ 
(187rV82(QM(z )-l)-39(QM(z)-l)^)/ ^M(z) using the spherical 
collapse model ( Brvan & Normanlll998i) for flat cosmologies. 
The mass of the halo is: 



where p,n - Pc ■ ^Atiz) is the mean density at the cluster 
redshift. We solved for the mass of t he cluster with the ex- 
pression for the shear yrir) derived by 'Bartelmann ("1996^ and 
Wright & Brainerd (2000) and minimized Equation[T4]with the 
MINUIT package. 

Table |2] shows the best-fit parameters and the mass values de- 
rived by model fitting. The NFW profile was used keeping 
both the concentration and mass as free parameters (marked as 
NFW), and by using the relation between Cyir and Mvir proposed 
bv lBullocketaPJMl (hereafter MNFW): c,,> = if:(^f, 

where M. = 1.5 x 10^^ /h Mq,K^9 and a = -0.13. 
These values were computed assuming spherical symmetry for 
the cluster halo. The effect of departures from spherical symme- 
try (e.g. triaxial halos) on the determination of the total clus- 
ter mass have been^tudied by several authors (e.g. Gavazil 
l2005l iDe Filippis et al.ll2005h . iOe Fifippis et all (l2005b showed 
that these efifects are negligible when the mass is computed at 
large distances from the cluster center, although they are impor- 
tant at small radii. The same authors tried to recover a three- 
dimensional reconstruction of Abell 611 through a combined 
analysis of X-ray and Sunyaev-Zel'dovich observations and con- 
cluded that the cluster was approximately spherical, supporting 
our symmetry assumption. 

In Figure [14] the results of fitting by a NFW model (dashed 
lines), a constrained (M)NFW {solid lines) model (iBullock et al.l 
2001) and a SIS profile {dotted lines) are displayed, for the KSB 
{left panels) and Shapelets {right panels) pipeline, respectively. 
The black lines represent the best fit to the unbinned data as func- 
tion of distance from the center of the cluster. Average values 
of tangential {black) and radial {red) components of the shear, 
computed in logarithmic scale bins, are also plotted. The latter 
components are expected to be zero in the absence of systemat- 
ics errors. 

In Fig.[T5]the confidence contours for the NFW profile are plot- 
ted in a plane r^ vs c: the levels show confidence at 68% and 
90%, starting with the innermost one. 

The goodness of each fit {y^) and its probability {Q) are also 
listed in Table |2] 



M,,,-,. = -7TAyirP,„r^,i^. 



(19) 
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KSB 


Shapelets 


Parameter 


MNFW 


NFW 


SIS 


MNFW 


NFW 


SIS 


M2oo(10") 
r2oo (kpc) 


1513+"'' 
4.51 


1545+3« 

T Q+5.6 




5 3+"->* 

-0.8,, 

1516!^^ 
4.50 


1570+'" 
3 T-^ 

'-1.3 




(Tci (km/s) 




778!55 






78 1!?^ 


red-x'' 


1.46 


2.02 


1.94 


1.62 


2.12 


2.31 


Q 


0.21 


0.11 


0.10 


0.17 


0.10 


0.06 



Table 2. Mass values computed for M200 by model fitting using a SIS profile, a NFW profile and a constrained NFW (MNFW), 
according with iBuUock et al.l (1200 ih . Best fit parameters and reduced listed for each fit performed using both KSB and 

Shapelets shear catalogs, together with the goodness of fit probability Q. 




Fig. 14. Results of fitting by a NFW model {dashed lines), a constrained (M)NFW {solid lines) model (iBuUock et al.ll200"l1) and a 
SIS profile {dotted lines), with the KSB {left panel) and Shapelets {right panel) pipeline. Average values of tangential {black) and 
radial {red) components of the shear, computed in logarithmic scale bins, are also plotted. Cluster masses were computed at r2oo, 
estimated to be ~ 350 arcsec. 



7. Discussion 

We have conducted a weak lensing analysis of Abell 611 with 
deep g-band images from the LBC. Due to the complexity of 
the LBC PSF, we decided to use both the KSB and Shapelets 
methods to measure galaxy shapes and extract the shear signal. 
KSB parametrizes the sources using their weighted quadrupole 
moments and is based on a simplified hypothesis of a nearly 
circular PSF. In contrast, Shapelets uses a decomposition of the 
images into Gaussian-weighted Hermite polynomials and does 
not make any assumption about the best PSF model. Due to the 
large collecting area of LBT and the wide field of the LBC, we 
were able to extract a high number density of ~ 25 background 
galaxies per arcmin^ over a wide field of 19 x 19 arcmin^ and 
this allowed us to perform an accurate weak lensing analysis. 
The two shear catalogs, derived by the KSB and Shapelets 
pipelines, were matched and common sources (with a number 
density of ~ 23 galaxies/arcmin^) were used to estimate mass 
measurements for Abell 611 by two different weak lensing 
techniques: aperture densitometry and a parametric model 
fitting. In both approaches we assumed that the BCG was the 
center of the mass distribution, which is supported by S-Maps 
(see Fig. [12] Section l6Tt . 

The projected mass values obtained by aperture densitometry 
within a radius of ~ 1500 kpc are: 7.7 ± 3.3 x IO'^^Mq and 



8.4 + 3.8 X W'^Mq, using KSB and Sha pelets, respectively . 
These estimations are model independent dCIowe et alJll998l) . 
but they are affected by large uncertainties. As discussed in 
Section|62] the contribution of uncorrelated structures along the 
line of sight can be source of noise for aperture measurements, 
although the y can be decreas ed by averaging the results for sev- 
eral clusters (lHoekstrall200II) or corrected by using photometric 
redshifts of the sources, if available. 

Table |2] shows the results of fitting the observed shear with a 
parametric model. We assumed both a SIS and a NFW mass 
density profile. We first fitted a NFW model leaving both c and 
rs as free parameters: as displayed in Fig. [15] the uncertainty 
on these two parameters provided by the fit is high and does 
not allow to put a strong constrain on the concentration. 
Nevertheless, the best-fit value (c ~ 4) is in good agreement 
with the one obtained when the Bullock et al. (2001) relation is 
adopted (c - 4.5). For the NFW fits the quoted masses are M21M, 
within the radius r2oo where the density is 200 times the critical 
density. The estimated value of r2oo for Abell 61 1 is ~ 1 .5 Mpc 
(5.8'). 

The weak lensing mass measurements obtained from both 
the KSB and Shapelets shear catalogs are in agreement, 
within uncertainties: by using a (M)NFW profile we obtain 
M200 ~ 5.3!j-^x 1O"*M0 and 5.3+JJ-^ x 10"*Mo, respectively. The 
smaller uncertainties of Shapelets results show that this method 
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Fig. 15. The confidence contours for NFW profile are plotted in the (rj,c) plane, for KSB (left panel) and Shapelets {right panel), 
respectively. The confidence levels are at 68% and 90%, starting with the innermost one. The triangle shows the best-fit position. 



provide a higher accuracy than KSB. Moreover, the goodness 
of the fits ix^) in Table 12 shows that both Shapelets and KSB 
provide a best-fit with higher probablity (Q) by using a NFW 
mass density profile than a SIS profile. 

Abell 611 was previously targeted for a weak lensing study 
by Dahle (2006), who used V and / observations from several 
facilities to target a large number clust ers (see [D ahle et al.] l2002l 
for more details). They used the KSB (iKaiser et al... 19951) shear 
estimator as described in lKaised(l2000l) to derive the shear signal 
from the images. The mass of the cluster was then derived by 
fitting the observed shear wit h a NFW mass density profile 
(H avarro. Frenk & Whitelll997l). Thev assumed a concentration 
parameter as predicted by Bullock et al.l (1200 ih and obtained a 
mass value of Mi8o = (5.21 + 3.47) X \Qi^%~^Mq within ri8o, the 
radius within which the density is 180 times the critical density. 
A more recent m ass es timate of Abell 611 is presen ted in 
[Pedersen & Dahle (l2007h . who used the data collected in iDahld 
(|2006) and obtained a value of M500 = 3.83 + 2.89 x IO'^/i-'Mo 
within rgoo, the radius within which the density is 500 times 
the critical density. In these papers the authors extrapolated 
the NFW profile out to rjoo because the data were insuffi- 
cient to extend this far in projection from the cluster center. For 
their work on Abell 6 1 1 r/^/rjoo - 0.59 (note rsoo - 0.66xr2oo). 



Our weak lensing estimates for the mas s of A b ell 61 1 are 
in agreement with the p revious results of iDahld (l2006h and 
iPedersen & Dahld (l2007h . but the depth and the larger area 
covered by LBC data allowed us to perform a more accurate 
analysis. 

A recent weak lensing ana l ysis of Abell 611 has r ecently 
been done by lOkabe et all (12009. submitted to PASJh using 
Subaru/Suprime-Cam observations in two filters (/' and V). The 
authors used the color (V - information to select the back- 
ground galaxies to use for their cluster lensing analysis, getting 
a galaxy number density of ~ 21 galaxies/arcmin-, with a SNR 
~ 10. By fitting the mass density profile of the cluster using a 
NFW model, they found a mass value Maoo 



and a c,,,. 



4 14+1.V3 
^•^^-1.21 



5.13![;i[;xlOi4M0 
(private communication), which are in 



iNewman et al.l (l2009l) over wide range of cluster-centric 
distance (from ~ 3 kpc to 3.25 Mpc) by combining weak, 
strong and kinematic analysis of the cluster, based on Subaru, 
HST and Keck data, respectively. They found a mass value 
M200 = 6.2+°-^ X IO'^Mq with c = 6.95 + 0.41 by using a NFW 
model fitting, in agreement with our results. We note that such a 
large value of c cannot be rejected by our NFW fits, due to the 
large uncertainties on NFW parameters (see Tab|2]and FigfTSll. 
Finally, we also compared the mass values obtained by our 
weak lensing analysis to X-ra y estimates of M200 available 
in literature. Schmidt & Allen I (2007) analyzed Chandra data 
of several clusters and modelled their total mass profile (dark 
plus luminous matter) using a NFW profile. They found for 



Abell 611 a scale radius rj - 0.32_y2o Mpc and a concen- 
tration parameter c - 5.39^[ jj, which provide a mass value 
A/200 ~ 8 X 10''* Mo at r2oo ~ 1.7 Mpc, in agreement, within 
the uncertainties, with the results obtained by us. A more 
recent X-ray analysis of Chandra observations of Abell 611 
has been performed by Donnarumma et al. (in preparation). 
They obtained a value of M200 - l-H ± 0.21 x IO'^Mq at 
r2oo ~ 1900 kpc, with r, = 407+if and c = 4.76+JJ-^^. Their 
projected mass at r2oo is 1.28 ± 0.24 x IO'^Mq. Such value is 
in agreement, within the statistical uncertainties, with the mass 
measured by aperture densitometry, but higher than the value 
estimated by the parametric model. Additional information on 
the mass will be derived from a strong lensing analysis of Abell 
611 (Donnarumma et al. in preparation). 

This work shows that LBT is a powerful instrument for weak 
lensing studies, but we want to stress that the data here analyzed 
did not allow us to use the full capabilities of the telescope. The 
presence of bright saturated stars hampered to use the whole 
field of the camera for the analysis of Abell 611. In addition, 
the analysis of weak lensing is expected to be improved by the 
usage of the Red Channel in LBC, which was not yet available 
during these observations. The present results are nevertheless 
important to demonstrate the capabilities of LBC for weak 
lensing; we therefore plan to extend to other clusters such 
analysis, now using the Red Channel. 



good agreement with our results. 
A new mass estimation of Abell 611 was performed by 
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